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ENSEMBLE DE TRAITEMENT THERMIQUE DE TISSUS BIOLOGIQUES ET 
PROCEDE DE MISE EN CEUVRE DE CET ENSEMBLE. 

^invention conceme le domaine des appareils destines aux therapies par 
5 hyperthermic locale. L'invention concerne egalement un precede de mise en ceuvre de 
tels appareils. 

Les therapies par hyperthermic locale consistent a chauffer, localement, une 
zone cible d'un tissu biologique. Lorsque ce type de therapie est utilise dans le cadre de 
la therapie genique, la chaleur peut etre, par exemple, utilisee pour son action sur un 
10 promoteur thermosensible. On peut egalement utiliser la chaleur pour necroser des 
tissus biologiques et proceder a 1* ablation de tumeurs. 

Aussi, les therapies par hyperthermic locale offrent-elles de nombreux 
avantages. Ces avantages sont aussi bien qualitatifs qu'economiques. Du point de vue 
qualitatif, elles ofifrent, par exemple, un fort potentiel pour le controle de traitements tels 
15 que les therapies geniques, le depot local de medicaments, l'ablation de tumeurs, etc. 
Du point de vue economique, elles sont compatibles avec un traitement ambulatoire des 
malades, elles permettent de reduire les durees dliospitalisation, etc. 

Dans les therapies par hyperthermic la chaleur peut etre, par exemple, apportee 
par un laser, des micro-ondes ou des ondes radiofrequences, des ultrasons focalises, etc. 
20 D'une maniere generate, les therapies par hyperthermic locale permettent des 
interventions medicales dont la nature invasive est reduite au minimum. Mais parmi les 
types d'energie precites, les ultrasons focalises sont particulierement interessants 
puisqu'ils permettent de chauffer la zone de focalisation, de maniere non invasive, 
profondement dans un corps biologique, sans chauffer de maniere significative les tissus 
25 voisins de la zone de focalisation. 

Dans tous les cas, la temperature de la zone cible et de son environnement 
immediat, pendant le traitement, doit etre controtee precis&nent et de maniere continue, 
alors que Tapport en energie est localise et rapide (de Pordre de quelques secondes). A 
cette fin, on peut installer des sondes de temperature dans la zone cible et son 
30 environnement immediat. Mais on peut egalement utiliser l'Imagerie par Resonance 
Magnetique (IRM). L'IRM permet, en effet, d^btenir une cartographie precise des 
distributions de temperatures ainsi que des informations anatomiques detaillees. L'IRM 
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permet en outre un controle non invasif de la temperature. 

On connait deja des dispositifs pour le controle de la temperature, pendant les 
traitements par des ultrasons focalises, bases sur la thermometrie par resonance 
magnetique. De tels dispositifs sont en particulier decrits dans les documents suivants: 

5 "Control system for an MRI compatible intracavitary ultrasons array for thermal 
treatment of prostate disease", Smith NB et al., Proceedings of the annual meeting of the 
International Society of Magnetic Resonance in Medicine, 1999, p. 672 et "Real time 
control of focused ultrasound heating based on rapid MR thermometry", Vimeux FC et 
al., Invest Radiol 1999, 34> p. 190-193. Dans les dispositifs decrits dans ces 

10 documents, le retrocontrole de la chaleur apportee par les ultrasons focalises, grace aux 
cartographies obtenues par IRM, est de type PID (acronyme de Pexpression angio- 
saxonne Proportional Integral and Derivative). De plus, avec ces dispositifs, le controle 
de la chaleur fournie au tissu, repose sur la prise en compte d'une temperature mesuree 
dans la zone de focalisation des ultrasons uniquement, ou correspondant a une moyenne 

15 obtenue a partir de la distribution spatiale de la temperature, dans la zone cartographies 
La figure 1 represente revolution temporelle de la temperature moyennee de la 
zone de focalisation, traitee grace au dispositif decrit dans le premier de ces documents. 
Sur cette figure, la temperature evolue jusqu'a un plateau correspondant a la 
temperature que Ton souhaite atteindre dans la zone de focalisation. On peut constater 

20 que la temperature souhaitee dans la zone de focalisation n'est atteinte qu'apres une 
duree de Tordre de 30 minutes. 

La figure 2 represente 1'evolution temporelle de la temperature moyennee dans la 
zone de focalisation, traitee grace au dispositif decrit dans le deuxieme des documents 
cites ci-dessus. On peut constater que la temperature souhaitee dans la zone de 

25 focalisation est atteinte en moins de 2 minutes. Par contre, on observe des variations de 
la temperature souhaitee, de plus ou moins 4° C. 

Un but de 1'invention est de proposer un ensemble pour le traitement thermique 
d'une zone cible d'un tissu biologique permettant a la fois d'obtenir rapidement la 
temperature souhaitee dans la zone cible et de maintenir et controler la temperature dans 

30 cette zone cible avec une precision accrue, par rapport a ce qui etait possible avec les 
techniques de Tart anterieur. 

Ce but est atteint, selon 1'invention, grace a un ensemble pour le traitement 
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thermique d'une zone cible d'un tissu biologique, comprenant 

- des moyens generateurs d'energie pour foumir localement de Fenergie dans la zone 
cible ; 

- des moyens pour mesurer et enregistrer la temperature dans la zone cible; 

5 - une unite de controle comprenant des moyens pour determiner, a partir de la 
temperature mesuree dans la zone cible, la quantite d'energie devant etre foumie a la 
zone cible, et des moyens pour commander aux moyens generateurs de Fenergie, de 
delivrer cette valeur de la puissance ; 

caracterise par le fait que 1' unite de controle comprend en outre des moyens de 

10 traitement numerique, point par point, de la distribution spatiale de la temperature dans 
la zone cible et son environnement, pour calculer des gradients de temperature. 

U ensemble de traitement thermique conforme a la presente invention, prend en 
compte la distribution spatiale reelle de la temperature dans la zone cible, mais aussi 
dans renvironnement de cette zone. C'est a dire qu'il prend en compte et traite, point 

15 par point, cette distribution spatiale. Contrairement aux ensembles de traitement 
thermique de Tart anterieur, la distribution spatiale de la temperature est utilisee pour en 
deduire des gradients de temperature et non pas de simples moyennes. Ceci permet 
d'estimer, avec plus de precision, quelle doit etre la quantite d'energie a deposer et done 
de parvenir plus rapidement a la temperature souhaitee et de maintenir la temperature du 

20 tissu biologique avec une plus grande stabilite. 

Avantageusement, Funite de controle de l'ensemble de traitement thermique 
selon T invention, comprend en outre des moyens pour estimer les pertes locales en 
energie thermique, a partir d'une estimation de la conduction thermique et de la 
distribution spatiale de la temperature dans la zone cible et son environnement. En effet, 

25 T information fournie par la valeur des gradients de temperature, ainsi que la prise en 
compte d'une estimation des pertes locales de chaleur, permettent non seulement de 
connaitre la fa<?on dont le tissu biologique traite a reagi a la chaleur dej& appliquee a 
celui-ci, mais permettent en outre, grace a la prevision sur la fa$on dont le tissu 
biologique va reagir a la chaleur. Ceci permet egalement de faire dvoluer la temperature 

30 du tissu traite thermiquement, plus rapidement vers la temperature souhaitee et de 
maintenir la temperature du tissu biologique avec une plus grande stabilite. 

Avantageusement, les moyens generateurs d'energie de Pensemble de traitement 
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thermique selon Finvention, emettent des ultrasons focalises. En effet, les ultrasons 
focalises permettent de fournir de la chaleur, dans une zone tres localisee, de maniere 
non invasive, meme si cette zone est situee profondement dans un corps humain ou 
animal. En outre la focalisation permet de ne pas chauffer de fa<?on significative, les 

5 tissus avoisinant la zone de tissu biologique traitee. 

Avantageusement, les moyens pour mesurer et enregistrer la distribution spatiale 
de la temperature de Fensemble de traitement thermique selon Finvention, comprennent 
un appareil d'imagerie par resonance magnetique. En effet, FIRM permet une mesure 
de la temperature, non invasive, precise, et bien resolue en de nombreux points de la 

10 zone cartographiee. Les donnees recueillies par IRM sont, en outre, aisement traitees 
numeriquement. 

Avantageusement, Tensemble de taitement thermique selon l'invention comprend 
des moyens devaluation de la distribution spatiale, dans la zone cible et son 
environnement, de l'energie fournie a la zone cible. 

15 Selon un autre aspect, Finvention est un precede de regulation d'un ensemble de 

traitement thermique d'une zone cible d ! un tissu biologique, comprenant 1'etape 
consistant a deposer localement de l'energie dans la zone cible, 
caracterise par le fait qu'il c -r- prend en outre les etapes consistant a 
- evaluer les gradients de temperature dans la zone cible et son environnement ; et 

20 - en deduire Tenergie a deposer dans la zone cible pour atteindre la temperature 
souhaitee. 

Avantageusement, ce precede comprend alors en outre Fetape consistant a 
estimer les pertes locales en energie, dans la zone cible et son environnement. 

Avantageusement, ce precede comprend en outre F6tape consistant a 6valuer la 
25 distribution spatiale, dans la zone cible et son environnement, de Fenergie fournie a la 
zone cible. 

D f autres aspects, buts et avantages de Finvention apparaitront plus clairement a 
la lecture de la description detaillee qui suit. 

Uinvention sera aussi mieux comprise a Taide des dessins joints sur lesquels: 
30 - la figure 1 represente revolution temporelle de la temperature d'une zone cible, 
lorsque celle-ci est traitee en utilisant un ensemble de traitement thermique de Tart 
anterieur ; 
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- la figure 2 represente revolution temporelle de la temperature cTune zone cible traitee 
par un autre ensemble de traitement thermique de Tart anterieur ; 

- la figure 3 represente schematiquement P ensemble de traitement thennique conforme 
al'invention ; 

5 - la figure 4 est un organigramme du procede de regulation conforme a la presente 
invention ; 

- la figure 5 represente revolution du Laplacien en fonction du temps, au cours du 
traitement thermique par 1'ensemble conforme a la presente invention, correspondant 
a une experience in vitro decrite ci-dessous ; 

10 - la figure 6 presente les resultats d'une serie d' experiences visant a etudier l'influence 

des erreurs d'estimation sur les pertes d'energie, representees par la diffusivite de la 

chaleur, et le coefficient d'absorption de Penergie ultrasonore ; 
_ la figure 7 represente revolution temporelle de la temperature maximum, mesuree au 

cours du traitement thermique d'un tissu biologique, par rensemble conforme a la 
15 presente invention, correspondant a la meme experience in vitro que celle de la 

figure 5 ; 

- la figure 8 represente revolution temporelle de l'amplitude du signal emis par le 
generateur de rensemble conforme a la presente invention, au cours du traitement 
thermique d*un tissu biologique, correspondant a meme experience in vitro que celle 

20 des figures 5 et 7 ; 

- la figure 9 represente 1'evolution temporelle de la temperature maximum, au cours du 
traitement thermique, d'un tissu biologique, par rensemble conforme k la presente 
invention, correspondant a une autre experience in vitro, avec trois paliers de 
temperature, et decrite ci-dessous ; et 

25 - la figure 10 represente Involution temporelle de la temperature maximum, mesuree 
au cours du traitement thermique par l'ensemble conforme a la presente invention, 
correspondant a une experience in vivo decrite ci-dessous. 

Un des modes de realisation de l'invention, est decrit ci-dessous de mantere detaillee. 
A titre d'exemple, ce mode de realisation de Tinvention correspond a un ensemble de 
30 traitement par hyperthermic locale par ultrasons focalises, controles par IRM. 
Comme represente sur la figure 3, un tel ensemble comprend : 
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des moyens generateurs d'energie 100 ; 

- des moyens de cartographie 200 ; 
une unite de controle 300 ; et 

- un porte-6chantillon 400 pour le tissu biologique 410a traiter. 

5 Dans le mode de realisation de F invention decrit ici, les moyens generateurs 

d'energie 100 sont composes d'un transducteur 110, d'un generateur de signal 
sinusoidal 120, d'un amplificateur 130 et d'un convertisseur 140, reliant le generateur 
de signal sinusoidal 120, a l'unite de controle 300. 

Le transducteur 110 opere a 1,45 MHz. Un tel type de transducteur 1 10 est, par 
10 exemple, commercialise par la societe Speciality Engineering Associates® (Soquel, 
Californie). Son diametre et sa distance focale sont de 38 mm et 25 mm respectivement 
Le generateur de signal sinusoidal 120 est, par exemple, du type FG110, 
commercialise par la societe Yokogawa® (Tokyo, Japon). 

L 5 amplificateur 130 est, par exemple, du type KMP 170F, commercialise par la 
15 societe Kalmus® (Bothell, Washington). Get amplifiacteur 130 a un gain en puissance 
de 58 dB. 

Le convertisseur 140 est, par exemple un convertisseur IEEE488 serie, 

commercialise par la societe I. O. Tech.® (Cleveland, Ohio). 

Les moyens de cartographie 200, permettent de mesurer et enregistrer la 
20 distribution spatiale de la temperature. lis comprennent, par exemple, un appareil IRM 

de type Bruker Biospec commercialise par la societe Bruker® (Ettlingen, Allemagne). 

Cet appareil utilise un aimant de 4,7 T qui est equipe d'un insert de 120 mm de 

diametre, qui genere des gradients du champ magnetique (la valeur maximum du 

gradient est de 0,193 T/m). 
25 L'unite de controle 300 comprend en particulier une station de travail 310 de 

type Alpha PW 500a MHz, commercialisee par la societe Digital®. 

L'unitg de controle 300 comporte egalement des moyens devaluation et de 

traitement numerique de la distribution spatiale de la temperature 320, des moyens de 

determination de la valeur de la puissance 330 devant etre fournie a la zone cible, des 
30 moyens d'estimation des pertes locales de Fenergie thermique 340 et des moyens de 

commande 350 des moyens generateurs d'energie 100. Les moyens de commande 350 

indiquent aux moyens generateurs d'energie 100 de delivrer la valeur de la puissance 
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foumie par les moyens de determination de la valeur de la puissance 330. 

Le porte-echantillon 400 comprend un support 420 de rat, en plexiglas®. Ce 
support 420 contient le transducteur 1 10 et une bobine de surface (non representee sur la 
figure 3). Un tel porte echantillon 400 a deja ete decrit dans les documents « Fast lipid 

5 supressed MR temperature mapping with echo-shifted gradient echo imaging and 
spectral-spatial excitation », de Zwart JA et al., 1999, Magn. Res.Med, 42, p.53-59 ; et 
« On the feasibility of MRI-guided focused ultrasound for local induction of gene 
expression », Madio DP et al., 1998, J. Magn. Res. Imaging!, 8, p. 101-104. Le support 
420 est place dans un tube de plexiglas® qui est partiellement rempli d'eau. Le 

10 transducteur 110 est positionne de maniere a ce que le point de focalisation 460 des 
ultrasons soit situe approximativement a 10mm en profondeur, dans le tissu biologique 
410. Pendant les mesures in vitro, une sonde 430 de temperature est inseree dans le tissu 
biologique 410 constitue d'un morceau de viande fraiche, de maniere a avoir une 
reference pour la temperature. Cette sonde 430 est par exemple un thermocouple du 

15 type Digi-Sence DualLog, commercialise par la societe Cole-Parmer Instrument Co.® 
(Vernon Hill, Illinois). 

La preparation des echantillons pour les experiences in vitro et in vivo est 
realisees de la maniere suivante. On prend des rats males de la race Wistar de 325 a 500 
g. Ceux-ci sont anesthesies en associant 1% en volume d'halothane avec un melange 

20 constitue de 7 volumes d'oxyde nitreux (N 2 0) pour 3 volumes d'oxygene, selon un 
protocole approuvS. Afin d'ameliorer la penetration du faisceau d'ultrasons dans le tissu 
biologique 410, la cuisse du rat (qui se trouve entre transducteur et point focal) est 
6pil6e, en utilisant un produit prevu a cet effet et connu de l'homme du metier. Pendant 
les mesures in vivo, la temperature endorectale des rats est enregistrSe. La temperature 

25 corporelle des rats est maintenue k 35°C, par immersion du corps des rats dans un bain 
dont la temperature est r^gulee a cette fin. Apres le traitement thermique les rats sont 
sacrifies. 

Les moyens de cartographie 200 sont utilises en mettant en oeuvre ime sequence 
Echo gradient, avec les parametres suivants : TR = 50 ms, TE = 15 ms, la taille de la 
30 matrice = 64x63, trois lignes de Tespace des k par TR, FOV = 64 mm x 63 mm, ou TR 
est le temps de repetition, TE est le temps d'echo et FOV est le champ de vue (FOV est 
Tacronyme de l'expression anglosaxone "Field of View"). Les donnees sont obtenues 
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par IRM, a partir d'une coupe de 2 mm d'epaisseur, perpendiculaire au transducteur a 
ultrasons et comprenant le point focal des ultrasons focalises. La resolution temporelle 
des cartographies obtenues par ERM est de 1,05 s. La resolution spatiale des 
cartographies obtenues par IRM est de 1x1x2 mm 3 . Les cartographies de la temperature 

5 mesuree par resonance magnetique sont obtenues a partir des mesures du deplacement 
de la frequence du proton de l'eau. Le choix de la resonance du proton de l'eau, pour 
faire ces mesures, repose sur le fait que la relation entre la frequence de resonance du 
proton de l'eau et la temperature est, en premiere approximation, independante de la 
composition du tissu biologique 410. Le deplacement de la frequence du proton de l'eau 

10 en fonction de la temperature est en outre lineaire et cette linearite n'est pas affectde par 
les modifications du tissu biologique 410 induites par la chaleur (Ishihara Y et al., 
Magn. Res. Med., 1995, 34, p. 814-823 ; Peters RD et al., Magn. Res. Med., 1998, 40, p. 
454-459). La dependance en temperature de la frequence de resonance du proton est de 
0,0094 ppm-K" 1 (« Fast magnetic-resonance temperature imaging » ;de Zwart JA et al, 

15 1996, J. Mag. Res. B, 112. p. 86-90 ; « Fast lipid supressed MR temperature mapping 
with echo-shifted gradient echo imaging and spectral-spatial excitation », de Zwart JA 
et al., 1999, Magn. Res. Med. t 42, p.53-59). 

Les signaux de resonance magnetique provenant des lipides constituent une 
source significative d'erreurs dans les cartographies de temperature calculees, puisque 

20 les frequences de resonance des protons des lipides ne dependent pas de la temperature. 
On supprime done les signaux de resonance magnetique provenant des lipides, en 
utilisant une excitation selective de l'eau, de la maniere decrite dans le document « Fast 
lipid supressed MR temperature mapping with echo-shifted gradient echo imaging and 
spectral-spatial excitation », de Zwart JA et al., 1999, Magn. Res.Med, 42, p.53-59. 

25 L'utilisation des moyens devaluation et de traitement numerique de la 

distribution spatiale de la temperature 320, a egalement deja ete decrit dans le document 
«Fast lipid supressed MR temperature mapping with echo-shifted gradient echo 
imaging and spectral-spatial excitation », de Zwart JA et al., 1999, Magn. Res. Med, 
p.53-59. 

30 Un exemple de mise en oeuvre du procede de regulation de Tensemble de 

traitement thermique, conforme a la presente invention est decrit de maniere detaillee 
ci-dessous. 
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Sur la figure 4, est represents schematiquement, un organigramme de ce mode 
particulier de mise en oeuvre du precede selon T invention. 

Conformement a cet exemple de mise en oeuvre, le precede selon P invention 
comporte : 

5 

- une etape 1 d'estimation du coefficient de diffusion de la chaleur (cti) et du 
coefficient d'absortion des ultrasons focalisSs (ot 2 ) dans le tissus biologique 410 ; 

- une etape 2 de definition par l'utilisateur, d'un profil de revolution temporelle 
de la temperature souhaitee ; 

10 - une etape 3 d'acquisition d'une image IRM ; 

- une etape 4 de calcul de la distribution spatiale de la phase au point focal 460 et 
dans son environnement ; 

- une etape 5 d'&ablissement de la distribution spatiale de la temperature au 
point focal 460 et son environnement ; 

15 - une etape 6 devaluation des gradients de temperature du point focal 460 et 

dans son environnement ; 

- une etape 7 de determination de la nouvelle puissance devant etre delivree par 
le generateur 120 ;et 

- une etape 8 de changement du taux d'energie deposee par le generateur 120. 

20 Les etapes 3 a 8 sont reprises en boucle, afin d'atteindre et de parcourir le profil 

de Involution temporelle de la temperature souhait6e, d^fini a l'etape 2. 

Pour un element donn£ du transducteur 110, la puissance electrique P(t) 
transmise a un echantillon est determinee par les moyens de determination de la valeur 
de la puissance 320. Sa valeur peut done etre modifiee en direct par Tunitd de controle 

25 300. Elle est obtenue sur la base de Tequation : 

Pit) = -^rj^ 1 ~ «.fcJ - Vr^o + a • [0(0- 7^(01+4 * A(0 J 

Cette equation est etablie en se fondant sur les considerations suivantes. 
Le point focal est d6fini par r= (0,0,0). devolution temporelle et la 
temperature maximum au point focal correspond alors a 7^(0= 7(0,0,0,0- Notons 
30 par 0(0 le profil predetermine de revolution temporelle souhaitee de la temperature 
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maximum T^Q). Comme indique ci-dessus, ce profil est defini avant le debut de 
chaque experience. Par exemple, ce profil ©(/) comporte une montee de 10° C, pendant 
100 s. Cette etape de montee suit revolution d'une moitie de la periode de la fonction 
cosinus. Eile est suivie par une periode ou la temperature est constante (10° C au-dessus 

5 de la valeur de depart), pendant 250 s. 

La temperature maximum 7^(0 = 7(0,0,0,0 P e **t etre controls seulement au 
point focal. En effet, la geometrie du transducteur 110 et la distribution spatiale de 
Findice de refraction dans le tissu biologique 410 determine le champ acoustique. En 
consequence, revolution de la temperature, ailleurs que dans la zone de focalisation, fait 

10 intervenir des fonctions dependantes de la coordonnee d'espace r et de la temperature 
T. Le champ de la puissance acoustique fij ), le tenseur de la diffusivite de la chaleur 
c^Qyt) dans le tissu biologique 410 et le coefficient d'absorption des ultrasons 
focalises a 2 (r , T) 9 sont alors lies par la relation : 

^^^(rjJ.Vn^o}^^^'^).?^). [la] 

15 

Lorsque la diffusivite est isotrope et varie lentement spacialement, l'equation la 
sesimplifie pour donner l'equation lb, ou a l (r 9 T) = ^Tr[a } (r,T)] est un champ scalaire 

d 1 d 1 & 

et V 2 est l'operateur Laplacien defini par —7 + —7 +-tj : 

3c dy azr 

= a, (r,T)V*T(r,t)+ a 2 (r,T)- p{r)- P(t) [lb] 

ot 

20 II faut noter que les fonctions o,(r,r) et o^if.T) ne sont pas precisement 

connues au debut du chauffage. 

Le champ de la puissance acoustique d'un element spherique de transducteur a 
ultrasons focalises, correspond approximativement a xme distribution gaussienne autour 
du point focal 460, avec un rayon d'attenuation a 6dB note R 0 . Le temps de diffusion 

25 specifique rest defini par r=-^-. Pour une longueur d'onde ultrasonore de 1 mm, 

1'ordre de grandeur de rest de l'ordre de 10 s. 
On se fixe les objectifs suivants : 
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1) la temperature du point focal 460 doit atteindre la temperature souhaitSe, aussi vite 
que possible (c'est-a-dire dans un temps du meme ordre de grandeur que r), ceci sans 
oscillation, ni depassement de la valeur souhait6e ; 

2) une fois cette valeur souhaitee atteinte, la temperature doit rester constante pendant 
5 une periode predefinie par l'utilisateur. 

L'integrale du profil 0(0, Tintegrale de revolution temporelle de la temperature 
maximum observee experimentalement T^t) et la difference entre ces deux integrates, 
sont respectivement : 

Q(,) = '$G(t)df [2a] 
o 

io ¥(0-/^(0* [2b] 

o 

a(0 = Q(0- *(>) = f r^O*)]* [2c] 

0 

En utilisant ces expressions, Tequation lb, donnant revolution de la temperature 
au point focal, peut etre exprimee en fonction de A(f): 

£^ = ^-a ] (T^).V^(t)-a 2 (T Bm )Pi0 [3] 
at at 

15 our = (0,0,0) est omis et p(0,0,0) = 1 . 

Dans Tequation 3, le parametre que Ton souhaite controier directement est la 
puissance des ultrasons focalises P(t). Notons qu'une equation differentielle d'ordre 2 et 
lineaire en A(f), peut etre avantageusement utilis6e par l'unite de controle 300, de 
maniere similaire a un systeme de controle PID. La raison en est que la solution pour 

20 A(/) d'une telle equation tend asymptotiquement vers zero, et qu f il en est de meme pour 
sa dSrivee premiere. Si la derivee premiere de A(f) est 6gale a zero, 7^(0 presente un 
recouvrement avec le profil predetermine de revolution temporelle de la temperature 
0(/). Ceci constitue l'id6e fondamentale du procede de controle mis en ceuvre par Vxm\\& 
de controle 300. Ainsi, nous pouvons reecrire liquation 3 sous forme d'une Equation 

25 differentielle lineaire de second ordre en A(t) de type : 
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d 2 A dA a 2 r41 
— r + a— +— -A = 0 15 J 

dt 2 dt 4 

De maniere a obtenir l'expression voulue de l'equation 4 a partir de l'equation 3, 
P(t) est reecrit avec l'expression suivante : 

TO = ^jff 4 " «.&.)• ^(0 + « • [©(0- r_Cr)] + 4 • [5] 

5 L'equation 5 correspond a l'equation centrale utilisee pour calculer en direct le 

niveau de puissance des ultrasons focalises. A partir de la solution de l'equation 
differentielle du second ordre, correspondant a l'equation 4, on peut voir que le 
parametre a est relie au temps de reponse caracteristique t r de la boucle de regulation, 
par l'expression a=2l t r . On suppose que dans les equations 4 et 5, toutes les fonctions 

10 utilisees pour calculer la puissance P(t) sont precisement connues. On peut egalement 
verifier, comine le montre Tequation 6 ci-dessous, que la temperature observee 
experimentalement T(t) 9 tend asymptotiquement vers le profit 0(0 : 

0(/)-r(O=f (0=[e(o)-r(o)].(i-^/)ex P (-|.r) ^ 

Comme deja note plus haut, dans une experience reelle, les coefficients 
15 d'absorption des ultrasons a> et le parametre de diffusion de la chaleur a/, ainsi que leur 
dependance en temperature, sont inconnus. Ces parametres et a 2 dependent de la 
composition du tissu biologique 410, de processus physiologiques, telle que la 
perfusion, et de changements irreversibles intervenant pendant la procedure de 
chauffage, par exemple dans les procedures d'ablation. Ainsi, un systeme de regulation 
20 doit il etre tolerant vis a vis des erreurs sur les parametres a, et a 2 . 

Seul le profil ©(/) et sa derivee sont precisement connues. 

Aussi, lorsque la puissance ultrasonore est calcutee en direct a partir de 
Tequation 5, deux difficultes apparaissent : 

1) 7^(0 et V 2 ^^) telles qu'elles sont obtenues a partir de la cartographie de la 
25 temperature issue de TIRM, sont aflfectees par le bruit. 

2) les valeurs de a, et a 2 et leur dependance en temperature ne sont pas precisement 
connues, de meme leur sensibilite a la necrose par chauffage (par exemple dans les 
procedures d'ablation) et les parametres physiologiques telle~que la perfusion, ne sont 
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pas connus avec precision. 

Toute erreur pouvant affecter a, et peut etre traitee comme une erreur de 
parametre dans une boucle de controle, selon un modele lineaire. Ainsi, des estimations 
des valeurs initiales deoj et seront choisies, puis utilisees pendant la procedure de 
5 chauffage, pour le calcul de la puissance ultrasonore, conformement a Tequation 5. 

Une analyse theorique de Peffet de rerreur sur les parametres a, et dans 
l'equation 5 met en evidence les effets suivants : 

1) Une estimation erronee du param&re diminue l'efficacite de la boucle de 
regulation ; il peut en effet s'ensuivre une possibility de depassement de la 
10 temperature souhaitee ou une valeur de la temperature determinee trop faible ; ceci a 

alors pour consequence une augmentation du temps de convergence. Quoiqu'il en 
soit, meme dans ces conditions, la temperature experimentale tend toujours 
asymptotiquement vers le profil predetermine de revolution temporelle de la 
temperature. 

15 2) Une estimation erronee du parametre a, conduit a un decalage constant, dans la 
zone ou le profil ©(/) est plat, des valeurs de la temperature, entre les valeurs de la 
temperature mesurees experimentalement et le profil 0(f). Ce decalage est 
proportionnel a la derivee premiere par rapport au temps, du Laplacien, multipliee 
par rerreur absolue sur a, , fois a' 2 . Afin d'estimer cet effet, la derivee du Laplacien 
20 a ete determinee en utilisant une regression lineaire de la courbe representee sur la 

figure 5, entre 150 et 250 s. Sa valeur est de Tordre de 0.01 K.mm'ls' 1 , ce qui 
conduit a un decalage en temperature approximativement egal a 0.1 °C. Ainsi, 
rerreur sur la temperature reelle ne devrait pas etre directement observ6e du feit de 
la limitation de la precision des mesures thermometriques par resonance 
25 magnetique, a cause du bruit. 

L'influence de rerreur sur les valeurs a x et a 2 sur l'efficacite du controle de la 
temperature par rensemble de traitement thermique conforme & 1'invention est 6tudi6e 
experimentalement, sur un morceau de viande fraiche, ci-dessous. 

Selon Tetape 1, les parametres a, et a 2 doivent etre estimes. Ceci est realise, 
30 dans un premier temps a partir d'une experience preliminaire avec une puissance 
d'ultrasons focalisds constante. Le parametre a x est calcule a partir de la derivee par 
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rapport au temps, de la temperature au point focal, divisee par la moyenne (sur cinq 
images IRM) du Laplacien. II est calcule immediatement apres ^extinction de 
remission des ultrasons focalises et est exprime en mm 2 /s. Le parametre a 2 (la vitesse 
de depdt de l'energie au point focal prenant en compte la distribution spatiale de la 

5 puissance des ultrasons focalises) est calcule a partir de la derivee, au temps initial, de la 
temperature du point focal, par rapport a la puissance des ultrasons focalises (lorsque les 
ultrasons focalises sont <§mis et que la diffusion est negligeable, voir ^equation 1). II est 
exprime en K^-^mV)" 2 . La precision estimee pour a, et est meilleure que 10%, 
comme on peut le deduire d'experiences repetees. Les valeurs numeriques obtenues de 

10 cette fa9on peuvent etre directement utilisees dans Tequation 5 pour calculer, a chaque 
fois qu'une nouvelle cartographie de la temperature est disponible, la valeur reelle de la 
puissance devant etre fournie par le generateur 120. 

Selon Tetape 2, un profil de revolution temporelle souhaitee de la temperature 
au point focal est defini avant le debut de chaque experience. Ce profil comporte une 

15 partie initiale croissante, correspondant a une moitie de la periode de la fonction 
cosinus, suivi d'une partie constante en temperature. La derivee premiere de la courbe 
correspondant a ce profil, est continue et peut etre calculee numeriquement par 
ordinateur. 

L'etape 3 est mise en oeuvre grace aux moyens de cartographie 200. 

20 Les etapes 4 a 6 sont mises en oeuvre grace aux moyens devaluation de 

traitement numerique de la distribution spatiale de la temperature 320. Les phases sont 
calculees a partir des signaux IRM, obtenus grace aux moyens de cartographie 200. Les 
changements de frequence de resonance de Teau sont calculees a partir de ces phases. 
Les changements de temperature sont calcules a partir de ces changements de 

25 frequence. 

La temperature maximum T^ft) et Tintegrale A(t) sont directement deduites des 
cartographie obtenus par IRM, respectivement au cours des etapes 5 et 6. L'integrale 
A(t) (Equation 2c) est calculee numeriquement a Taide de la station de travail 310. 

Au cours de l'etape 6, les pertes locales de Tenergie thermique sont reevaluees a 
30 Taide des moyens d'estimation des pertes locales de Tenergie thermique 340. 

Selon Tetape 7 du precede conforme a la presente invention, Foperateur 
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Laplacien V 2 7*(r,/) est applique au traitement des cartographies de la temperature 
obtenues par IRM. La valeur de ce Laplacien au point focal est calcule, en utilisant la 
methode des elements fins, en combinaison avec un filtre temporel de reduction du 
bruit. Ce filtre utilise un rapport de ponderation binomial egal a 1 :4 :6 :4 :1, sur cinq 
5 images. Cette etape 7 est mise en oeuvre avec les moyens de determination de la valeur 
de la puissance 330. 

Finalement, la puissance est calcutee grace a l'equation 5 et la valeur calculee est 
envoyee au convertisseur 140 (etape 8). L'etape 8 est mise en oeuvre avec les moyens de 
commande 350 des moyens generateurs d'energie 100. 

10 La duree totale de calcul pour traiter chaque cartographie de distribution spatiale 

de la temperature, c'est a dire chaque cycle de Tensemble des etapes 3 a 8 decrites ci- 
dessus, est inferieure a 250ms. 

Afin d'analyser la tolerance vis-a-vis des erreurs dans les estimations initiates de 
Oj et plusieurs procedures de chauffage ont ete mises en oeuvre avec les parametres 

15 «j etcr 2 variant sur une large gamme. Cette gamme va de 0 a 300 % d'une valeur 
predeterminee, pour a, , et de 40 a 250 %, pour . Cette valeur predeterminee est celle 
obtenue a partir des mesures preliminaires exposees ci-dessus. Des temps tfattente de 30 
minutes ont ete introduits entre les experiences successives dans le but d'atteindre les 
lignes de base en temperature, spatialement uniformes dans Techantillon et identiques 

20 pour chaque experience. Neuf resultats representatifs sont reportes sur la figure 6. Ces 
resultats montrent que le syst&ne a une large tolerance vis-a-vis des erreurs sur les 
estimations des valeurs initiales de a, et Oj. Quoi qu'il en soit, on peut noter que la 
boucle de controle devient instable seulement lorsque est fortement surestim^. Cette 
instability est exacerbee lorsque a, est fortement surestimd et a 2 est sous-estim6e. Cet 

25 effet peut etre attribu^ au bruit experimental des mesures de la temperature par IRM. 
Lorsqu'un calcul affecte par le bruit, conduit a une surestimation de la valeur du 
Laplacien (les derivees secondes etant sensibles au bruit), la puissance des ultrasons 
focalises appliques augmente dans une proportion egale a aj.£fa 2 , ou s est la 
surestimation du Laplacien. Une augmentation de la puissance des ultrasons focalises 

30 conduit a une forte augmentation du Laplacien dans le tissu biologique 410 et de ce fait, 
a une nouvelle augmentation de la puissance des ultrasons focalises. Cette reaction 
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positive cessera apres un temps approximativement egal a 2/a, grace a une reaction 
negative dans la boucle de controle. Ceci explique la periodicite des instabilites dans ce 
cas extreme. 

La force de la reaction negative sur la boucle de regulation correspond au 
5 parametre a. En effet, comme nous Favons vu plus haut le parametre a, est egal a deux 
fois Tinverse du temps de reponse caracteristique t r (a=QJ f r ). La valeur de ce parametre 
est indiquee pour chaque exemple de la figure 6. En general, les valeurs de a de 0,1 s" 1 & 
0,2 s' 1 sont suffisantes pour atteindre des temps de montee en temperature similaires a 
ceux du profil ©(*)> meme lorsque des valeurs fortement erronees de aj et o^sont 

10 utilisees. Cest seulement dans le cas extreme ou aj est fortement sous-estime et o& est 
surestime (voir figure 6, en bas a droite), que a doit etre augmentee jusqu'a la valeur de 
0,40 s" 1 pour obtenir un recouvrement avec le profil predetermine de revolution 
temporelle de la temperature. La valeur de a optimale, experimentalement trouvee (dans 
tous les cas, sauf les cas d'erreurs extremes sur et Oj), est de 0,2 s" 1 , le temps de 

15 reponse de la boucle de regulation correspondante etant de 10 s. Lorsque la force de la 
reaction negative est augmentee, une correction plus rapide des erreurs sur les 
parametres initiaux est obtenue, mais Tamplitude de la puissance des ultrasons focalises 
et des fluctuations en temperature autour de la valeur predeterminee est aussi 
augmentee. 

20 La mise en oeuvre et les performances de F ensemble de traitement selon 

T invention decrit ci-dessus, sont iliustrees ci-dessous a Faide de deux exemples. 

Exemple 1 : Utilisation de Fensemble de traitement thermique, conforme a la 
presente invention, dans le cadre de mesures in vitro 

Selon cet exemple, un protocole de mon.ee en temperature de 10°C est mis en 

25 oeuvre sur un echantillon de viande fraiche. La temperature initiale est egale a 15°C. 
Avec ce protocole, aucune modification irreversible du tissu biologique 410, resultant 
de ce profil devolution temporelle de la ter :rature, n'est attendue. La figure 7 
represente revolution de la temperature maximum en fonction du temps. Dans la partie 
plate de cette courbe, Felevation moyenne de la temperature est de 9.97°C, avec un 

30 ecart type egal a 0.1 9°C. Cet ecart type doit etre compare a celui, egal a 0.1 8°C, qui est 
obtenu par les mesures de temperature effectuees sans chauffage par ultrasons focalises 
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(c'est a dire ce qui correspond a la ligne de base du bruit dansrles mesures de 
temperature). Sur la figure 5, est reporte le Laplacien calcule en direct. L'attenuation 
observee sur la partie constante en temperature correspond a la diminution des gradients 
en temperature autour du point focal. L'amplitude de la puissance appliquee en direct 

5 est reportee, quant a elle, sur la figure 8. A cause du bruit de mesure, la valeur calcul6e 
du Laplacien et l'amplitude de la puissance fournie par le generateur, presentent une 
fluctuation de Pordre de 10% , approximativement Ceci n'a qu'un faible effet sur la 
temperature resultante, car la frequence de fluctuation (c'est a dire l'inverse de la 
resolution temporelle de la cartographie par resonance magnetique) est beaucoup plus 

10 grande que l'inverse du temps de reponse (r) specifique au chauffage du tissu biologique 
410. 

La figure 9 montre la stabilite en temperature obtenue avec un profil a trois 
paliers (15, 25, 30 °C). L'ecart type est egal a 0.35°C, 0.36°C et 0.40°C, respectivement 
pour des montees en temperature a 15°C, 25°C et 30°C. Les resultats representes sur la 

15 figure 9 confirment la grande stabilite en temperature, sur une large gamme de montee 
en temperature, du systeme de regulation de 1'ensemble de traitement thermique 
conforme a la presente invention. 

Exemple 2 : Utilisation de 1'ensemble de traitement thermique, conforme a la 
presente invention, dans le cadre de mesures in vivo. 

20 En adoptant une procedure analogue a celle mise en ceuvre dans le cas de 

T exemple 1, des experiences ont ete menees in-vivo, sur une cuisse de rat. Les resultats 
correspondants sont montres sur la figure 10. La resolution temporelle est de 0.5s. La 
temperature moyenne, entre 90 et 120 s apres le debut de Pexp^rience, est de 54.9°C (la 
valeur du profil a atteindre est 6gale a 55°C) avec un ecart type de 0.33°C. Les figures 

25 7, 9 et 10 montrent que Ton peut controler la temperature avec une precision proche de 
celle donnee par les mesures de temperature effectuees in vitro ou in vivo. 

II a ete decrit ci-dessus un mode de realisation de Tinvention correspondant k un 
ensemble de traitement par hyperthermic locale par ultrasons focalises, controles par 
IRM, mais T invention couvre une gamme plus large d' ensembles de traitement 

30 thermique. Aussi, on comprendra que Tinvention peut etre generalisee aux cas ou la 
chaleur est, par exemple, apportee par un laser, des micro-ondes ou des ondes 
radiofrequence, des ultrasons focalises, etc. On comprendra egalement que d'autres 
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moyens de mesure de la temperature, peuvent etre utilises dans Fensemble de traitement 
thermique conforme a la presente invention, a la place de FIRM. 

De meme revaluation et la traitement numerique de la distribution spatiale de la 
temperature ont ete decrits ci-dessus comme etant effectues a Taide du Laplacien. 
5 D'autres moyens d'effectuer cette evaluation peuvent etre utilises sans sortir du cadre de 
Tinvention ; 
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REVINDICATIONS 

1. Ensemble pour le traitement thermique d'une zone cible d'un tissu biologique 
(4 1 0)fcomprenant 

5 - des moyens gen^rateurs d'energie (100) pour fournir localement de Tenergie dans la 
zone cible ; 

- des moyens (200)pour mesurer et enregistrer la temperature dans la zone cible; 

- une unite de controle (300) comprenant des moyens (330) pour determiner, a partir de 
la temperature mesuree dans la zone cible, la quantite d'energie devant etre fournie a la 

10 zone cible, et des moyens pour commander (350) aux moyens generateurs de Tenergie 
(100), de delivrer cette valeur de la puissance ; 

caracterise par le fait que l'unite de controle (300) comprend en outre des moyens de 
traitement numerique (320), point par point, de la distribution spatiale de la temperature 
dans la zone cible et son environnement, pour calculer des gradients de temperature. 
15 2. Ensemble de traitement thermique selon la revendication 1, caracterise par le 

fait que l'unite de controle (300) comprend en outre des moyens (340) pour estimer les 
pertes locales en energie thermique, a partir d'une estimation de la conduction 
thermique et de la distribution spatiale de la temperature dans la zone cible et son 
environnement. 

20 3. Ensemble de traitement thermique selon Tune des revendications precedentes, 

caracterise par le fait que les moyens generateurs d'energie (100) emettent des ultrasons 
focalises. 

4. Ensemble de traitement thermique selon Tune des revendications precedentes 
caracterise par le fait que les moyens (200) pour mesurer et enregistrer la distribution 

25 spatiale de la temperature comprennent un appareil d'imagerie par resonance 
magnetique. 

5. Ensemble de traitement thermique selon Tune des revendications precedentes, 
caracterise en ce qu'il comprend des moyens devaluation de la distribution spatiale, 
dans la zone cible et son environnement, de l'energie fournie a la zone cible. 
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